Effects of Quantum Error Correction on Entanglement Sudden Death 
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We investigate the effects of error correction on non-local quantum coherence as a function of time, 
extending the study by Sainz and Bjork [1[ . We consider error correction of amplitude damping, pure 
phase damping and combinations of amplitude and phase damping as they affect both fidelity and 
quantum entanglement. Initial two-qubit entanglement is encoded in arbitrary real superpositions 
of both <3>-type and ^P-type Bell states. Our main focus is on the possibility of delay or prevention 
of ESD (early stage decoherence, or entanglement sudden death). We obtain the onset times for 
ESD as a function of the state-superposition mixing angle. Error correction affects entanglement 
and fidelity differently, and we exhibit initial entangled states for which error correction increases 
fidelity but decreases entanglement, and vice versa. 

PACS numbers: 03.67.Pp, 03.65.Ud, 03.67.Bg, 03.65.Yz 



I. INTRODUCTION 

Entanglement is a condition shared by vectors in two 
or more separate vector spaces. As is well known quan- 
tum entanglement can provide an important resource for 
quantum computation and communication because it is a 
non-local effect, i.e., it can exist among multi-party sys- 
tems completely disconnected from each other (see, e.g., 
Preskill jij and Nielsen and Chuang Q). The quantum 
entanglement in a multi-party system is a form of joint 
coherence that is subject to change whenever the system 
interacts with an external environment, typically becom- 
ing smaller. Such coherence-diminishing interactions are 
almost always present, and since they can reduce or elim- 
inate the usefulness of entanglement their effects are im- 
portant to understand. 

Environmental interactions can cause loss of coherence 
in single systems as well as joint systems, but entangle- 
ment is a more fragile coherence in some ways Q. A 
particular example is provided by what is labelled ESD 
(early stage decoherence or entanglement sudden death) . 
This can take place in common physical contexts [f| even 
under relatively benign dissipative interactions with the 
environment. 

As the name implies, ESD is more destructive than 
single-system decoherence, and not in the sense that it is 
governed by a larger rate constant (which is usually not 
the case), but because ESD completely eliminates entan- 
glement in a finite time (see manifestations in various 
model contexts: [6 10] and the review [111]). This can 
happen while the local coherences of the same parties re- 
main continuously non-zero under the same environmen- 
tal interaction. ESD is not well understood. Its onset 
is not predicted in any known way for a general entan- 
gled mixed state (however, see recent results [Hj]). This 
is true even for the simplest case, qubit-qubit entangle- 
ment, which can be present in the smallest possible (4 x 
4) joint density matrix. 



Quantum error correction 13] is a known weapon 
against single-system decoherence, but tests of its efficacy 
against multi-party decoherence and specifically against 
dissipation of quantum entanglement are rare [if, [lij ]. 
Here we report results of calculations that expand the do- 
main for which results are available. We focus on the ef- 
fects of quantum error correction on ESD under different 
noisy- environment interactions, using error correction al- 
gorithms that are convenient for this limited purpose. In 
practice a noisy channel is a dynamical process and its 
physical features are usually not fully understood or pre- 
dictable, and so are potentially a complication for error 
correction. Additionally, the multi-qubit logic states that 
are required to introduce needed redundancy multiply 
the overall resource requirements. 

We particularly extend the quantum error correction 
study initiated by Sainz and Bjork [l[ in which they 
showed that error correction can actually promote ESD 
instead of preventing or delaying it. We present compar- 
isons of time-dependent results for uncorrected and error- 
corrected state evolution, focusing on two-qubit entan- 
glement with concurrence as measure of entanglement. 
A comparison is made of concurrence with fidelity, also 
frequently used to describe multiparty coherence. 

The method that we use lets the error-channel in ques- 
tion affect our system for an amount of time before cor- 
rection. That is, we assume that the qubit in question 
represents a stored resource between protocol-controlled 
gate operations and it is during this "quiet" interval that 
most environment-induced errors occur. After this in- 
terval we isolate our system and subject it to error cor- 
rection. The later the isolation, the higher the chance of 
error occurrence. The error correction methods we adopt 
encode the qubits that are to be protected into multi-bit 
logical codewords and they can perfectly correct the er- 
rors that effect only one of the bits in these codewords. 
The encoding operation is local. Namely, given a pair 
of entangled and physically separated qubits, we apply 
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encoding onto each of them at their own stations. We 
will assume that this operation does not alter the entan- 
glement between the qubits as we will not try to extract 
information on the states of the qubits during the encod- 
ing. Finally, we will assume that the apparatus used in 
the error-correction is itself error-free. 

It is possible to use codes that correct multi-bit er- 
rors. However, one has to increase the number of bits 
in the codewords in order to be able to correct higher 
order errors. This in turn may increase the probability 
of uncorrectable errors depending on the probability of 
a single-bit error. For instance, it has been shown that 
Shor's 9-bit code [15| can be used to correct amplitude 



damping errors that effect up to two bits in the code- 
words [16]. Fig. [T] shows a comparison between this code 
and the 4-bit code we use. It is clear that unless the 
single-bit error probability, p, is smaller than about 0.1, 
the success probability of the 4-bit code is higher. Since 
we will be dealing with the whole range of error proba- 
bilities, < p < 1, the single-bit error correction is more 
appropriate in our case. 




FIG. 1: The probability of occurrence of a correctable error 
vs. p, the probability of an error at a given bit. The effective- 
ness of the 9-bit code decreases substantially with respect to 
the 4-bit code for larger p. 



II. AMPLITUDE DAMPING 

Amplitude damping here refers to the results generated 
by the interaction between a two-level quantum system 
(qubit) and a weakly coupled environment that is mod- 
eled as consisting of a continuum of harmonic oscillators. 
The boson modes of the electromagnetic vacuum are one 
example, but the model is practically general - any very 
weakly excited environmental interaction will behave as 
a harmonic oscillator. The damping effect is to drive 
the qubit to just one of its two possible states, thus re- 
moving any possibility of coherent superposition of the 
states. The Hamiltonian representing the interaction can 
be written, 



(1) 



where a± are the raising and lowering Pauli operators 
for the atom and a k , a\. are the annihilation and cre- 
ation operators for the k th mode of the environment. The 
strength of interaction between the atom and each mode 
is determined by the parameter g k . Under the Born- 
Markov approximation 17|, the master equation for the 
time evolution of the single-qubit density matrix is given 
in the interaction picture by, 



P = r ad (2cr -pa j, 



-P - 



(2) 



where T ad is the rate of damping which depends on the 
spectral properties of the environment and the interac- 
tion parameter. The Kraus operators representing this 
Markovian evolution are, 



E = 





V 1 ~ Pad 
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(3) 



where p ad = 1 
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evolution implied by Eq. 
the Kraus operators: 



2 can then be obtained using 



(4) 



where the subscript e indicates error-channel evolution. 
This Kraus representation ensures that the AD channel 
has no effect on the |0) state while it changes |1) into 
|0) with probability p a d and leaves it unchanged with 
probability 1 — p a d- This can be translated into a state 
amplitude-changing map for the joint atom-boson mode 
state: 



|0)|?mc) 
\l)\vac) 



\vac) 



VI -Pad \l)\vac)+y/p^\0)\<i>), (5) 



where |0) is the one-boson state created in the process. 

The error correction process works in four steps: First, 
the state is encoded using codewords that consists of mul- 
tiple bits. Second, the encoded state is exposed to the 
error channel. Third, the system is isolated and a syn- 
drome measurement is performed in order to find out the 
type of error that occurred. Finally, conditional on the 
result of the syndrome measurement, a recovery opera- 
tion is performed on the system in order to (partially) re- 
store the original state. Thus, given the initial encoded 
state po, and its error-modified form p e , the recovered 
state, pR is given by, 



p R = Y,RkM k p E MlRl, 



(6) 



where M k is the syndrome measurement operator that 
detects the k th type error and R k is the corresponding 
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recovery operator. We describe the method more fully in 
Sec. EHl at the end. 

Now we address non-local coherence under error- 
channel evolution and correction. Following [1], we apply 
this procedure to a pair of entangled (but possibly re- 
motely located) qubits, both interacting only with their 
local environments, taken to be identical and exposing 
them to the same type of amplitude-damping error. First 
we test the entanglement of a generalized Bell-type state 
|$), which is a superposition of the usual |$ ± ) states, 
known to be ESD-susceptible, to represent the initial 
atomic state: 



|$), a = 7r/4 



|$) = cos<S|$ + ) +sin<5|$ _ ; 
= cosajll) + sina|00), 



(7) 



where a + 5 = it /A. Since the |11) and 1 00) states fit the 
error correction steps better than the \^} states, we will 
call a the mixing angle for the superposition. 

We then calculate the fidelity and degree of entangle- 
ment of the corrected state and compare each with the 
values for the uncorrected one. We measure entangle- 
ment with concurrence and under the environment's 
amplitude-damping influence we find for concurrence and 
fidelity: 



C#(p e ) = max 0,2(1 -p od )| cos a| 



sina| - | cosa\p ad ) , 



(8) 



= 1 - 2p ad cos 2 a + p 2 ad cos 2 a. (9) 

There is no general analytical formula for C$(/jr) and 
F$(Pr), the concurrence and fidelity of the state after 
the error correction. However, it is possible to find such 
formulas with a small p ad approximation [l[ . 

Fig. [2]shows the results for maintenance of concurrence 
and fidelity for two values of the initial superposition 
mixing angle a. Time is represented parametrically by 
Pad = 1 — e~ r - 4Dt , ranging from to 1 (t=0 to t=oo). The 
figures should be interpreted as follows: Our qubits start 
from a pure state and evolve under the error channel. The 
dashed lines represent this evolution (we refer to them as 
"uncorrected"). For each time point, the corresponding 
point on a solid line represents the state that is reached 
from the dashed line after error correction. 

Several features are worth note in Fig. [2] For example, 
fidelity is confirmed as a bad substitute for concurrence, 
i.e., it is not a reliable measure of entanglement. Fi- 
delity is finite for all time, whereas entanglement vanishes 
at a finite time in all cases except the uncorrected pure 
Bell state (a — 7r/4). The remarkable feature first noted 
by Sainz and Bjork [1] is shown in the top left panel, 
where correction of the amplitude damping effect actu- 
ally initiates ESD rather than preventing it. However, 
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FIG. 2: Concurrence (top) and fidelity (bottom) vs. p for 
|$) under amplitude damping, showing the different effects of 
error correction depending on a — n/A (left) and a = 7r/12 
(right). The dashed line represents p e , the error-affected but 
"uncorrected" state, and the solid line represents pa which is 
the mapping of p £ under the error correction operation, thus 
referred to as the "corrected" state. 



one can also see that a departure from pure Bell status 
(a = 7r/12) leads to more rapid ESD, but the error cor- 
rection in that case is somewhat effective, not preventing 
ESD but at least delaying it. We can also see that the 
time of error correction must be chosen carefully. On the 
top left panel, it can be seen that around jad = 0.2 the 
error correction maps the error-affected state into one 
with a higher concurrence. However, if the error correc- 
tion is instead applied at jad = 0.4 the state is mapped 
into one with a lower concurrence. In fact, it is evident 
from the same figure that there is no use in applying er- 
ror correction once time reaches around jad — 0.3. In 
order to increase the entanglement, correction must be 
applied before this point. 

For the other Bell state we write its generalized super- 
position as 



I*) 



cos^|* + ) + sin 5|* 
cosa|10) + sina|01) 



(10) 



where again a + 5 = tt/A. We also calculate the fidelity 
and degree of entanglement of the corrected states and 
compare them in Fig. [3] with those of the uncorrected 
ones. Under the environment's influence we find: 



(AD) 
Fx (AD) 



2| sinacosa|(l 
1 -Pad, 



Pad) 



and 



(11) 
(12) 



so that fidelity is independent of the mixing angle, and 
entanglement of the state remains finite for all time 
(Pad < 1) in the absence of testing for error. 

However, Fig. [3] shows that when error correction is 
instituted (and amplitude errors are surely injected by 
the damping), one consequence is improvement in con- 
currence (top panels) for a very short period, but quickly 
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FIG. 3: Concurrence (top) and fidelity (bottom) vs. p for 
the |$) superposition under amplitude damping, showing the 
effects of error correction in its dependence on two superpo- 
sition angles, a — n/A (left) and a — n/12 (right). Error 
correction is of little consequence for fidelity. It provides pos- 
itive assistance to entanglement, but only for a small fraction 
of the amplitude damping halflife, after which it quickly in- 
duces ESD. 

such improvement disappears and a strong ESD event is 
induced by the error correction. As in the <£> state case, 
the response of fidelity to error correction again has little 
relation to the response of entanglement. 

III. PHASE DAMPING (PD) 

As a second error channel, we examine pure phase 
damping (PD), which in some situations can be more 
rapid than amplitude damping. It can be interpreted as 
arising from rapid stochastic shifts of the qubit transi- 
tion fre quen cy, which has a clear classical analog (see 
e.g., dHl), and can be easily extended beyond the 
Markov limit [2l| . Using bosonic reservoir language, in 
pure phase damping a boson is created or annihilated 
in the environment while the energy of the atom is un- 
changed; thus it can also be interpreted as arising from 
weak scattering processes. The interaction part of the 
Hamiltonian engages the qubit 's inversion, and is given 
by, 

H i D = °* Y^ hkak + h *k a i)- ( l3 ) 

k 

The master equation for phase damping in the inter- 
action picture can be written as, 

P = ^pd(<JzP<J z - p), (14) 

where Tpu is the phase-damping parameter. Only the 
off-diagonal elements are affected: 

Pn = Pn(0) and p 2 2 = P22(0) 

P12 = Pi2{V>)e~ 2VpDt and p 21 = p* 12 , (15) 



and the Kraus operators for this evolution are, 

*-(i^)'*-0!X)- <16 > 

where p p d = 1 — eT 2TpDt . 

Once again, this code can fully correct the first order 
errors, where one bit is flipped. In this case, the con- 
currences and fidelities for the corrected and uncorrected 
states are equal to each other: 

C*(p e ) = C^{PD) 

= 2| sinacosa|(l — jpo) (17) 

F*(p e ) = F^(PD) 

= 1 - 2j PD sin 2 a cos 2 a. (18) 

Note that uncorrected concurrence is above zero for all 
times (p p d < 1), only vanishing at p p d = l(t = oo). How- 
ever, the correction algorithm makes a difference. Fig. SI 
which serves for both |<£>) and |^), shows that error cor- 
rection (solid line) counteracts phase damping and main- 
tains the entanglement of a single Bell state (a = ir/A for 
either $ or "J/) above the uncorrected value for all time, 
but when applied to Bell superposition states (a ^ 7r/4), 
can give rise to ESD, as in the upper right panel. Sim- 
ilarly, fidelity of a single Bell state is well protected for 
arbitrarily long times, but is very poorly protected for 
a Bell superposition even for short times, as seen in the 
bottom right panel. 



$),a = tt/4 |$>,a = 7r/12 
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FIG. 4: Concurrence (top) and fidelity (bottom) are plotted 
vs. Ppd under phase damping for a — it/4: (left) and a — n/12 
(right). As before, the effect of elementary error correction is 
shown in blue. The same figures apply to both |3>) and 



IV. RELATIVE EFFECT (AD AND PD) 

Fig. [S] shows the effect of the error correction on the 
concurrence and fidelity more clearly. Here we introduce 
the differences between corrected and uncorrected quan- 
tities for a direct indication of correction effectiveness. 

Ac = Cv(p R ) - C 9 (ft e ) 
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(19) 



$) and |*) 



The interesting point is that there are regions for both 
values of superposition angle a where Ap is positive while 
Ac is negative. In other words it is possible that the 
error correction increases the fidelity while it decreases 
the entanglement. We will see in the next section that 
the opposite can also be true. 

From the data displayed above we can easily extract 
the absolute differences that are found between corrected 
and uncorrected records. These differences are shown in 
Figs. |5]and [5] for amplitude damping and phase damping 
respectively. Here we are defining differences as 



A c = C^( PR ) -C9(p e ) 

Aj? = F*(p R )-F*(p £ ). 



(20) 



The interesting point is that there are regions for both 
values of a where Aj? is positive while Ac is negative. 
In other words it is possible that the error correction 
increases the fidelity while it decreases the entanglement. 

Fig. [6] shows the effect of error correction on the con- 
currence and fidelity of the |$) state for two different a 
values. The right panel has regions where Ac is positive 
while Aj? is negative. In other words, error correction 
may lead to an increase in entanglement and a decrease 
in fidelity. 



V. COMBINED NOISE (AD AND PD) 

A remaining task is to examine combined AD and PD 
errors, since in general we would expect to have both of 
these error mechanisms present together to some extent 
in a physical system. The interaction part of the Hamil- 
tonian in the combined noise case is 




Pad 



FIG. 5: The differences between the concurrences and fideli- 
ties with or without correction for |3>) and j^) under AD. 
The results are plotted for two different values of a = n/4 
and q = 7r/12. 
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FIG. 6: The differences between the concurrences and fideli- 
ties with or without correction under the PD. Both |"1>) and 
j^) have the same dynamics under PD. 



jjAD+PD _ Y AD (2(7-p(T + - <7+<T-p - p<7 + <7-) 

+ T PD (a z pa z - p), (21) 



where Tad and Tpp, are the damping parameters for 
the AD and PD respectively when they are applied sep- 
arately. Clearly, the relation of the two damping rates 
may be important, and we will designate it by the ratio: 

K = T PD /T AD . (22) 

The solution for qubit evolution is given by, 



Pll = f-p22(0)e- 2r ^ 4 

PX2 = pi2{0)e- TAD - 2TpDt 

P21 = P*12 

P22 = P22{Q)e- 2TADt . 



(23) 



A two-operator Kraus sum for this type of combined error 
does not exist. Instead we can find a three-operator sum 
with the elements, 



1 

-Pad)(l -Ppd) 

1 




- Pad)Ppd 



(24) 



where p a d and p p d have already been defined. In order to 
make the AD and PD parts of the sum more transparent, 
we apply a unitary transformation to these elements to 
turn them into, 
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FIG. 7: Concurrence and fidelity are plotted vs. p a d for [$) 
under the combined noise for a — n/4 (left) and a = n/12 
(right) using k = 1 and re = 10 (ft = Fp^/Tad). F° r re = 10 
entanglement of both corrected and uncorrected states vanish 
before p ad — 0.4 whereas fidelity always stays above zero. 



correspond to the errors that cannot be corrected and 
if the syndrome measurement points to one of them we 
project our original qubit onto (|0)(0| + |l)(l|)/2. The 
results are given in Fig. [7] and Fig. [5J 

The concurrences and fidelities for the uncorrected 
states are given by, 

C$(p £ ) = max 0,2(1 -p Qfi )| cos a| 

xHsina|(l - p pd ) - |cosa|p afi ) , (28) 

F<s>(p e ) = 1 - tyad cos 2 a + p 2 ad cos 2 a 

-2p pd (l - p ad ) sin 2 a cos 2 a (29) 
C*(p £ ) = 2| sinacosa|(l - p a d){l - P P d) (30) 
F^(p £ ) = l-p ad -2p pd (l - p ad ) sin 2 a cos 2 a. (31) 

VI. ESD ONSET TIMES 



E = 
E' x = 

eL = 




(25) 



where, 



x = 1 + yj (1 -pad)(l -Ppd), 

V = \J (1 -Pad)Ppd- 



(26) 



E x and E 2 are the error operators for the PD and AD 
respectively. E represents the no-error case for small 
p ad and p pd - 

We will use the [5,1] code proposed by Laflamme et at 
[lj| for the error correction (see Sec. VII). As in [l|, we 
will ignore the distortion caused by the second diagonal 
element of E and assume that the no-error subspace is 
spanned by |0l) and The one- bit PD (AD) error 

subspace is spanned by the vectors each of which results 
from the application of £1(62) on the relevant qubit of 
\0l) and and identity operators on the rest, where, 



£1 



£2 



1 

-1 

1 





(27) 



In the end we get 22 orthogonal vectors, each correspond- 
ing to a different type of one-bit error (or a no-error) 
originated in |0l) or |1l)- The remaining 10 vectors that 
are necessary to span the 32-dimensional codeword space 



The factors that govern or control ESD are still gener- 
ally unknown, particularly in the case of initially mixed 
state evolution. One feature that would be valuable to 
have under control is the onset time for ESD. To our 
knowledge, no attempt to examine the onset time sys- 
tematically has been made. A recent exception is found 
in 



12] 



We now let T be the time when C reaches zero, i.e., 
T is the ESD onset time. Then from (j8j we obtain the 
ESD onset value of p a d for |<fr) as 



Pad(T) = I tana] 



(32) 



However, some results of this kind can be obtained 
from an extension of our examinations of concurrence for 



= tt/4 




k=1 , corrected 

K=1, uncorrected 

k=10, corrected 

k=10, uncorrected 
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FIG. 8: Concurrence (top) and fidelity (bottom) are plotted 
vs. p a d for \\V) under the combined noise for a — n/4 (left) 
and a = n/12 (right) using two different k values, k = 1 and 
ft = 10. Increasing k decreases the ESD onset time. Also, 
as for the AD and PD cases, error-correction leads to smaller 
ESD onset times. The features seen in Fig. 5 are observed 
here as well. 
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both uncorrected and corrected evolution. We have com- 
puted the evolution of concurrence, under the correction 
algorithms used above, for a wide range of a values. For 
essentially every a there is found a finite onset time for 
ESD. In Fig. |H] we show the way the ESD onset time, 
measured by the appropriate p, depends on a. 




FIG. 9: ESD onset time, measured by the relevant p, is plotted 
vertically vs. the Bell superposition parameter a for AD (left) 
and PD (right). Uncorrected does not show ESD under 
either AD or PD, and uncorrected |$) does not show ESD 
under PD. As a result, these cases are not plotted. 

The dashed line in the left panel of Fig. |H] shows that 
uncorrected concurrence is preserved for long times if the 
initial state was nearly a single Bell state (a w tt/4). The 
solid line shows that application of the error correction 
algorithm rapidly induces ESD for all values of a. In the 
right panel the solid line shows that phase damping is 
much more effectively forestalled by the correction algo- 
rithm used. Only values of a substantially different from 
7r/4 lead to a short ESD time under phase damping. 

It is evident that | combination does not show ESD 
under combined AD and PD effect. For |$), the ESD 
onset time is determined by, 



Pad(T) 

l-Ppd{T) 



tana 



(33) 



This should be compared with the ESD onset for AD 
which is given by, 



Pad{T) = | tan a | . 
Since for any given time, 



Pad 



1 



> 



P P d 



Pad, 



(34) 



(35) 



Eq. 38 is satisfied at an earlier time than Eq. 39 (except 
for a = 0). Thus, T decreases with the introduction of 
the PD to an AD-only system. The uncorrected states 



- k=1 , corrected 

- k=1 , uncorrected 

- k=10, corrected 

- k=10 uncorrected 




FIG. 10: ESD onset time for measured by p a d is plotted 
vertically vs. the Bell superposition parameter a for the com- 
bined noise using two different k values, k = 1 and k = 10. 
ESD onset time decreases with increasing n. For k = 10, 
uncorrected state has a larger ESD time than the corrected 
state for the whole range of a. Corrected ys>) gives the same 
results. Uncorrected Y&) does not show ESD. 



do not give ESD in a PD-only system. The introduction 
of even a slight amount of AD, however, leads to a finite 
Pad/(1 — Ppd) which means ESD will set in after a finite 
amount of time when AD is introduced into a PD-only 
system. 



VII. SKETCH OF METHOD 

The main challenge in the error correction process is 
finding a suitable code. The [4, 1] code given by Leung et 
al. [2 2| | is an example of such a code, and we have used it 
for amplitude damping. It encodes the two qubit states 
into two logical states that each live in a 4-dimensional 
Hilbert space. Each qubit is encoded using four bits: 



|o) L = (loooo) + |ini))/A 

\l) L = (|0011> + |1100»/V2, 



(36) 



and the scheme can fully correct only the errors that 
come from a single AD error which works like a one-way 
bit-flip, annihilating |0) while changing |1) into |0). The 
AD channel is then applied separately onto the four code 
bits which in itself gives an error operation consisting of 
16 operators (2 operators for each qubit). 

For example, suppose that we have a single-qubit state 
given by a\0) +b\l) which is encoded as o|0)l + 6|1)l and 
exposed to the AD error. Assume that, during the error 
process, an error occurs only at the first qubit in the code. 
Following the notation used by Sainz and Bjork [l|, the 
states that are reached from the codewords under this 



error are: 



|no) 



|0111> 
|0100), 



(37) 
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where the first index 1 shows the number of the qubit 
that is affected by the error and the second index or 1 
shows which codeword includes the affected qubit. Note 
that the two states are orthonormal. This is an essential 
property which is necessary to have, to be able to dis- 
tinguish the source of the error. It is also essential that 
the states for different types of errors are orthonormal 
as well. This ensures that one can distinguish between 
different error types upon a syndrome measurement and 
thus can recover the initial state correctly. 

Returning to our example, the syndrome measurement 
operator for this error type is given by, 

Mi = |0111) (0111 1 + |0100) (0100|, (38) 

and the corresponding recovery operator is given by, 

i?i = |0) l (0111| + |1)l(0100|, (39) 

so it can be shown that, 

R^hia^L + b\l) L ) = a\0) L + b\l) L . (40) 

Thus the original state is retrieved. The other one-bit 
errors and the no-error case are recovered similarly. It 
is shown in 1] that the vectors that are reached from 
the codewords under these errors are all orthonormal, so 
these errors can be distinguished from each other and 
corrected. 

The main problem with this procedure is that the cor- 
rectable errors constitute only a fraction of all the errors. 
For instance, in this case, there are 5 types of correctable 
errors but in total we have 16 error types. Some of these 
errors are not detected by this procedure at all and the 
ones that are detected get mixed up with one of the five 
errors that can be corrected. For instance the error type 
where the first, third and fourth qubits are damped lead 
to the vector |0100) when applied to \0)l- However, as 
we have seen above, the same vector is also reached from 
\1) l upon application of the error type where only the 
first qubit is damped. Thus, our procedure treats this 
error as if it stemmed from and recovers it as such 
while it should have actually given |0)l. The probabil- 
ity of having these high-order errors is small if time (or 
equivalently the damping parameter p a d) is small. Thus, 
the correction mechanism gives good results for small 
times. However, this weakness becomes more visible as 
time increases and the correction eventually breaks down, 
meaning that leaving the state uncorrected gives better 
longer-term entanglement and fidelity values. 

As for the PD case, a unitary recombination of Eq and 
Ei in Eq. [16] gives a new and equivalent set of operation 
elements for phase damping [see Theorem 8.2 of 0|], 

E'o = VP ( I I ) . E[ = ^T—p (\\), (41) 



where 

/3=i±4^. (42) 

One can see that the effect of the error operator E 1 above 
is to change the state |+) = (|0) + |l))/v2 into the state 
I - ) = (|0) — |l))/v2 an d vice versa. There is a three-bit 
error correcting code for the type of errors that flip the 
two states into each other in which our qubit states 
are encoded as 

|o)i = |-—>, 

= { + ++)■ (43) 

The states reached by a single-bit error from these code- 
words are orthonormal to each other and to the code- 
words. Thus, the code can accurately detect and recover 
from single-bit errors. However, it is not reliable for cor- 
recting higher order errors. For instance, an error that 
affects the first two qubits of |0)l is detected to be an er- 
ror that affects the third qubit of 1 1) x, and thus the state 
is recovered as |1) instead of |0). 

In the case of combined AD-PD damping we have used 
the [5,1] code proposed by Laflamme et al. [l9( for the 
error correction. The codes that we have used for AD 
and PD were specifically designed to cope with one-bit 
errors of their respective error type. This code, on the 
other hand, can correct general one-bit errors. In our 
case, it can be used to correct the errors where only one 
of the qubits in the codewords is effected by either AD 
or PD. Errors that effect more than one qubit cannot be 
corrected with this code, as it was the case for the codes 
given in previous sections. The logical bits of the code 
are given by, 

|o) L = Qooooo) + |11100) — |10011) — jOllll) 

+ |11010) + |00110) + |01001) + |10101)V2\/2, 
\l) L = (- |00011) + |11111) - 1 10000) + |01100) 

+ |11001) - |00101) - |01010) + |10110))/2V2. 

(44) 

There are 10 types of error that this code can correct: 
Half of them coming from the application of a single AD 
error and the other half from the application of a sin- 
gle PD error on each qubit in the logical bits. One can 
show that the states reached from the codewords through 
these errors are orthonormal to each other and to the 
codewords. As a result, given that one of these errors oc- 
curred, a syndrome measurement can determine exactly 
which has occurred and the initial state can be restored 
through the application of the proper recovery opera- 
tion. However, as in the [4,1] code, higher order errors 
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get mixed up with the correctable errors and this leads 
to a decrease in the efficiency of the code as the error 
probability increases. 

VIII. SUMMARY AND CONCLUSIONS 
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We have examined the dynamics of two-qubit systems 
under amplitude and phase damping, where the noises 
are both applied separately and as a combination, with 
and without the assistance of error-correction procedures. 
The examples are not complicated, and the error correc- 
tion algorithms are generically specific but not sophisti- 
cated. However, we believe they are general enough that 
conclusions can be drawn about preservation of prepared 
quantum coherence in arbitrary superpositions of generic 
|$) and \^f) Bell states. Coherence can be measured in 
several ways, and we employ both fidelity and concur- 
rence. 

It is well-known that concurrence fl8| is a completely 
reliable measure (both necessary and sufficient) to evalu- 
ate the degree of entanglement that may be present in a 
two-qubit state. Fidelity is examined because it is some- 
time used as a substitute in place of concurrence, and 
implications about entanglement are sometimes drawn 
from it. However, it is not a measure of entanglement, 
and not a reliable substitute. As we have seen, it suffers 
its own difficulties under error correction, not correlated 
in an obvious way with those for entanglement. 

The main conclusions are straightforward: As Sainz 
and Bjork [l| first discovered for amplitude damping, as 
we have confirmed in Fig. [31 an error-correction pro- 
cess can itself give rise to ESD, i.e., kill entanglement 
completely. This clearly establishes that error correction 
does not always preserve entanglement. Additionally, as 
in the combined AD and PD case, it can decrease the 
ESD onset time significantly, depending on the strengths 
of the AD and PD damping parameters. 

We do see that error correction can act positively with 
good success in some cases. Unfortunately, a systematic 
discrimination of cases is still lacking. This is illustrated 
for pure phase damping in Fig. 21 where concurrence is 
strongly protected for all p a d values in the case a = tt/4, 
while it is not for a = n/12, for which it goes into ESD 
after 2-3 half-lives of the dissipation. Second, a finite fi- 
delity does not necessarily indicate a finite entanglement, 
which is shown by amplitude damping for both a values 
tested. Overall, it seems safe only to say that error cor- 
rection may affect entanglement and fidelity in opposite 
ways, increasing one while decreasing the other. 
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